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Abstract 

We present a re-analysis of the cosmic confusion hypothesis, elucidating the degree 
to which "confusion" can be expected to hold in a class of flat, adiabatic models. 
This allows us to devise a simple and accurate fitting function for the height of the 
first peak in the radiation power spectrum in a wide range of inflationary inspired 
models. The range preferred by current data is given, together with a discussion of 
the impact of measurements of the peak height on constraining models of structure 
formation. 
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I. INTRODUCTION 



The picture of structure formation through gravitational instabihty now appears to be well 
established. Initial fluctuations grow through gravitational instability to form the structures which 
we observe today, and leave their imprint at redshift z ~ 1000 in the cosmic microwave background 
(CMB) radiation anisotropy. Within this framework any theory which purports to explain the 
large-scale structures we observe today, must simultaneously fit the increasing amount of data on 
CMB anisotropies. The combination of these two constraints is especially powerful, since they probe 
a large lever arm both in scale ik < 10~^h~^Mpc to k > l/i~^Mpc) and time {z > 1000 to 2; ~ 0). 

Calculations of both CMB anisotropy spectra and matter power spectra are now well developed 
for inflationary models, with the former approaching 1% precision [0,0. (In Fig. |^ we show the 
anisotropy power spectrum for the "standard" CDM model, see 0-^] for a discussion of the physics 
behind these anisotropies). This kind of precision in theoretical predictions is necessary to inter- 
pret the results envisioned from a future satellite CMB mission, which could accurately measure 
the power spectrum over a broad range of scales, allowing one to measure most of the standard 
cosmological parameters with unprecedented accuracy However such a measurement is still 

(optimistically) several years away, and in the meantime the data on CMB anisotropies continues 
to be amassed at a steady rate. For the interpretation of current data, such a high level of precision 
is not necessary, and semi-analytic techniques exist which can give spectra accurate to ~ 10% 0. 

An important constraint on structure formation models can be obtained without even this 
amount of effort however. It has been pointed out that for measurements of the spectrum 

constrained to scales of ~ 0°5 or above, the anisotropy power spectra exhibit a degeneracy in 
parameter space, which can be both a hindrance and a boon. This degeneracy has been dubbed 
"cosmic confusion" |^ as it provides a limit to how well cosmological parameters can be extracted 
from measures of the CMB power spectra on scales larger than 0°5. 

Notice the promise this degeneracy has in simplifying the interpretation of degree scale CMB 
data. As noted by Dodelson & Kosowsky 0, if this degeneracy holds it is possible to characterize 
the constraints on a wide range of models in terms of two parameters: an overall amplitude and the 
value of the degenerate combination of parameters u (see §|l|for details. A part of this approximate 
degeneracy was implicitly used in ^ to constrain the spectral slope and ionization history from 
the peak height). The amplitude of the potential fluctuations, measured on the largest scales by 
COBE, has become the preferred method for normalizing theories of large-scale structure |l^ . 
(We show an example of this in Fig. |l|, where the dimensionless amplitude of the matter power 
spectrum at horizon crossing, 6h, is shown vs. fio = 1 ~ and spectral tilt n.) The value of u 
will constrain other parameters on which the predictions of the theory depend. Furthermore, to 
the extent that cosmic confusion holds, generating spectra for a wide range of parameters can be 
reduced to simple modifications of a standard spectrum (however with the advent of semi-analytic 
techniques for quickly generating spectra over the entire angular range, this may be of limited 
utility). An obvious application of this is a fitting formula for the height of the peak in the power 
spectra at £ ~ 220. The current data provide limits on the height of the peak (see e.g. 0), which 
are becoming stronger. A quick method of estimating the height of the peak in the CMB spectrum 
can allow one to estimate the constraints from degree-scale measurements quickly and efficiently 
over a broad range of parameter space. 

The outline of this paper is as follows. In §^ we discuss how well cosmic confusion works, and 
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slightly modify the original formulation §|IT| discusses the special case of late reionization in 
processing the spectra. In § [rV| we use a compilation of CMB data to provide a limit on the peak 
height and discuss implications for models of structure formation. §|V| presents the conclusions. 



II. COSMIC CONFUSION 

In this section we would like to make some comments about, and minor improvements upon, 
the degeneracy of CMB spectra with respect to certain parameters which has been dubbed "cosmic 
confusion" . Specifically, it is claimed that for a range of parameters the CMB anisotropy spectra 
are degenerate in 

n = n- 0.28 ln(l + 0.8r) - 0.52 (nl^^h - O.s) - 0.000364^^ (1) 

where n is the (scalar) spectral index, r is the ratio of tensors to scalars (see below), h is the Hubble 
constant in units of 100km s~^ Mpc~^ and zr is the reionization redshift (see §|IITD. The degeneracy 
is supposed to hold up to the first peak in the power spectrum, i ~ 220, though some spectra can 
be degenerate over an even larger range of i 

To start let us take as our guiding principle the following conjecture: 

• For reasonable input parameters all CMB anisotropy spectra with the same peak-to-plateau 
ratio will be approximately degenerate. 

The interpretation of this statement clearly depends upon what range of cosmological models 
are adopted, over what range of scales the degeneracy is required to hold (not beyond the first 
peak), and the precise meaning of "approximately". We will discuss all of these issues below. 

Inflation in its most generic form predicts a fiat universe, so in this paper we will restrict 
ourselves to models with vanishing spatial curvature. Several questions relating to open inflation 
are still unresolved, for example the predictions for super-curvature modes, spectral tilt and tensor 
modesQ While for any given model the calculations of the power spectrum are straightforward 
[]T2|,3 , we feel that these uncertainties make fitting to spectra premature for these models. 

The simplest implementation of cosmic confusion is in terms of an apparent spectral index, 
which we shall call u (this is related to the of [§] by ~ n — 1). As our base model we shall 
adopt 'standard' CDM (h = 0.5, ^Iq = 1, f^A = 0, n = 1, = 0) for which we define u = 0. This 
spectrum is shown in Fig. |^ and tabulated in Table |. The spectrum for a model whose parameters 
combine to give u is then related to that of standard CDM by the simple formula 



The contribution to the large-angle CMB anisotropy from tensor perturations will probably be small in 
open models, since a combination of the COBE and cluster abundance normalizations predicts the spectral 
index n to be very close to (or greater than) 1. 
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^K^) = [y^J D,{u = 0) (2) 
where Di = £{i + l)Ce for notational convenience. We have chosen to pivot around £ = 10 as this 



is approximately the pivot point of the COBE data [0. Alternatively one could use cr(10°) 0, 
which turns out to be almost the same for the range of spectra under consideration. 

To understand the expression for u in terms of the parameters of the model, we need only consider 
the effect of changing each parameter on the peak height of the model. We will now consider each 
of these effects in turn, concentrating first on those parameters for which the conjecture works well 
and then turning to those for which the deviations are larger. 

As discussed in changing the spectral index, n, of the primordial spectrum (oc k^) is equivalent 
to multiplying the anisotropy spectrum by as in Eq. ^. For n > 0.8 the worst deviations induced 
by this approximation occur around i = 30-50, and are ^ 10%. Hence tilting the model away from 
scale invariance changes z/ by n — 1. 

Associated with tilt is the possibility that some of the CMB anisotropy comes from long- 
wavelength, inflation-produced gravitational waves (tensor fluctuations). The fraction, r, of tensors 
is usually quoted in terms of the contribution to the quadrupole: r = C2 /C*2 • Unfortunately 
due to the decaying potentials in a A model the scalar quadrupole is very sensitive to [0,|1^] 



for a fixed "initial" power spectrum. The tensor quadrupole is much less sensitive |T6| and hence 
the ratio has a strong (artificial) A dependence]^ jl^. We choose instead to write r = 1.4C{q /C{q , 
where the prefactor has been chosen to make the definitions agree for Qq = n = 1. The i = 10 mode 



is not as strongly affected by the decaying potentials |27] making this a more robust measure of 



the underlying tensor/scalar ratio. Since the tensor spectrum damps rapidly on scales smaller than 



the horizon at last scattering |T6| , p!9| the effect of introducing a tensor component is to increase the 
plateau by (1 -|- 0.76r) while leaving the peak unchanged. [The factor 0.76 is roughly / for 
a scale-invariant spectrum with fig = 1-] An equivalent reduction in the peak height relative to the 
plateau can be accomplished by reducing the spectral index (from £ = 10 to £ ~ 220) by, 

z/ 9 - ln(l + 0.76r)/ ln(220/10) - -0.32 ln(l + 0.76r) (3) 

In the semi-analytic, two-fluid models of CMB anisotropics pO|-|2^J^J^, the height of the first 
peak depends primarily on the combinations fi^/i^ (which sets the speed of sound in the coupled 
baryon-photon fluid) and the redshift of matter-radiation equality. The latter scales with Qoh'^f^, 
where = 1.68p^/prad = 3.36/(7* is the fraction of the radiation made up of photons, scaled to 
be 1 in the standard model with 3 massless neutrino species. In models in which more relativistic 
species are included, such as decaying neutrino models p3|, can be reduced leading to a larger 
peak height, just as if Qq were lowered. As a warning we point out that this 'degeneracy' is not exact, 
due to the effect of Q\ at low redshift (which is included implicitly in the formulation of the semi- 
analytic models [^]). In Fig. ^ we show spectra for two models with the same (reasonable) values 



^This dependence was originally neglected in |jT8| thus overestimating the damping effect of tensors for 
the tilted models. 
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of VLfiK^f^ and ^sh^, which differ even up to the first peak by 5-10%. Specifically the models have 
Qoh^ = 0.1 and Qeh^ = 0.0125, the first with ((]o,^^) = (0.4,0.25) and the second with (1,0.1). As 
one can see from the figure, the exact height and position of the first peak are complicated functions 
of the cosmological parameters! The reason for the breaking of the degeneracy in this particular 
case is that the large-angle temperature fluctuations (from decaying potentials at A domination) 
and the projection of physical scales to i both depend on Qa, not just (1 — VLi^h?. 

The original statement of cosmic confusion was meant to hold for models with Vtsh"^ = 0.0125, 
meaning that models with equal Vt^h? would be degenerate, and furthermore that they would 
coincide with models with some n ^ 1. For models with scale invariant initial perturbations, 
VLq = 1, Qsh'^ = 0.0125 and h = 0.3 to 0.75, we find that the height of the peak (relative to -Dio) is 
very well fit by 



D 
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We note that ipeak ranges from 260 to 200 over this range of h, since a changing sound speed causes 
the sound horizon at last scattering to subtend a varying angular scale This small movement 
of the peak should not matter if we redefine "cosmic confusion" to hold when averaged over a 
reasonably broad window in i, as is the case in most experiments to date. Under the "confusion" 
assumption then we predict that 

1/ 9 -0.371n(2/i) (5) 

which is roughly equivalent to the original statement given in 0. Both this approximation, and 
the original statement P|, give ^ 10% deviations in Dg for 10 < £ < 250 and 0.35 < h < 0.75. [The 
percentage deviations for £ < 10 can be larger than 10%, but here cosmic variance is also larger.] 

Similarly we find that the variation with fio = 1 — at fixed fi^/i^ and h works at the ~ 10% 
level in power for i > 10 and Qq > 0.4, with 

u 3 -0.161n(fio) - 0.15 ln(/^) (6) 



This is not exactly what we would have expected based on a degeneracy in ^Vt^f^h: the scaling with 
VLq is slightly weaker, and that with weaker still. [This can be traced to the same non- degeneracy 
mentioned above and illustrated in Fig. ^.] 

Note that at present VlBh'^ as determined by Big Bang Nucleosynthesis (BBN) is uncertain to a 
factor of ~ 2 (c.f. the Hubble constant!) ||2^. Thus we are compelled to study the variation of the 



peak height with the sound speed, or ilsh , which was not included in the original statement of 
cosmic confusion [^]. For fig = 1 and 0.01 < Qb < 0.10 the height of the peak is relatively well fit 
by an exponential in Qb for h = 0.5 and less well fit by an exponential for h = 0.8. Unfortunately 
the slopes of the fits for these two cases (which have different floh'^) disagree: 



In ( -^p^ ~ ig^Bh^ + 1.388 for h = 0.5 (7) 
~ 23^5/1^ + 0.824 for h = 0.8 (8) 
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As a compromise then we take the average of the two coefficients to arrive at our dependence of v 
on Q^B^ 



z/ 9 6.8(^^/1^ - 0.0125) (9) 



which holds relatively well over the range preferred by BBN: 0.01 < < 0.02 

In summary then we can define a "spectral tilt" v through 

V = n-\- 0.32 ln(l + 0.76r) + 6.8(^5/^2 - 0.0125) 
-0.371n(2/i) -0.161n(no) -0.151n(/^) 



(10) 



For a large range of Jltot = 1 CDM models, we find with this v that "cosmic confusion" holds at the 
~ 10% level in power (5% in temperature) up to £ ~ 200 for those models with fio > 0.4, h > 0.4 
and 0.01 < flsh^ < 0.02 (we discuss r > in §|TTp. This deviation usually occurs at £ ~ 30-50 
and some of it can be attributed to the approximation for tilt of the primordial spectrum that we 
are using (Eq. For low-i7o or h models the shape of the rise into the peak and the position of 
the peak are sufficiently different that a tilted model is not a good fit (though our tit is a good 
approximation to the height of the peak even for Qq ~ 0.3). On the other hand for some models 
in the range mentioned the degeneracy is good to 1% in power for 10 < £ < 200. The logarithmic 
dependence of u on Vt^h'^f^ differs from the linear dependence of h on the same quantity found by 
[^. Both approximations give similar "worst ffi" Di, but the logarithm is a better ffi to the height 
of the peak (-D220 relative to Diq) over the range of models considered; it works to better than 5%. 

With the rapid progress being made in measurements of the CMB anisotropy, it is likely that 
we will need to do better than 10% in the very near future. A more complicated fit or a better 
treatment of spectral tilt could alleviate matters slightly. However, bearing in mind that the position 

1 /2 

of the peak also changes with the input parameters and that the degeneracies (e.g. VLq h) built into 
"confusion" themselves only hold at the 5-10% level, one is lead to conclude that "one number" 
summaries of the CMB data are becoming a thing of the past. Of course, Eq. |10| is still useful as 
a quick-and-dirty method of estimating the height of the peak (to 5%) in the power spectrum for 
quite a wide range of parameters. This can be useful for narrowing the large parameter space down 
to a smaller region which can be searched more carefully (as one often uses linear theory estimates 
in large-scale structure work). If the experiments are chosen to probe scales for which "confusion" 
works well, then constraints on v will still encode much of the information that the CMB has for 
large-scale structure. 

As more experiments probe scales "beyond the ffist peak" (£ > 200), the utility of cosmic 
confusion as a summary of all CMB constraints will decrease. In these cases one must search 
a multi-dimensional parameter space and compute the full CMB spectrum either by numerical 
evolution of the coupled Boltzmann equations or fast semi-analytic methods which can be accurate 
to < 10% (e.g. 0). Moving beyond this degeneracy will allow us to obtain more information about 
the cosmological parameters from study of the detailed structure of the power spectrum. 
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III. REIONIZATION 



In the original statement of cosmic confusion |p, there was a term for the reionization redshift. 
Reionization in the adiabatic models of interest is likely to occur very rapidly |^ and at relatively 
low redshift, zr. Since zr the universe has probably been fully ionized. We show in Table |T| 
a "reasonable" estimate of the reionization redshift for some models by way of example. These 
determinations rely on several assumptions and so should only be taken as illustrative. The results 
are also very sensitive to the value of Vtsh? and chosen. 

In such a scenario, the amplitude of the fluctuations on small scales is reduced by e"^"^ ||2^j26|j2^ , 
where 



T = 0.035 -^hxe 



(for Qq + = 1) is the optical depth to Thomson scattering from z = to zr (also shown in 
Table |l[). Note that r depends not only on zr but also on Qb and h, which was neglected in 
the original treatment of cosmic confusion ^j. (For the accuracy to which the spectrum depends 
on and scales with r see f^.) In addition new fluctuations are generated on larger angular scales 



for T > 0.1 [22,2s]. Reionized spectra are not well described by a simple tilting of the sCDM 
spectrum. An approximation which has the correct asymptotic forms is to multiply the spectrum 
by exp[—2T{zi)], where t{z') is the optical depth from z = to redshift z = z' and zi defines a 
mapping from angles to redshift (e.g. the redshift at which the horizon subtends an angle 
While this approximation has the right asymptotic form, it fails to be a good approximation to a 
reionized spectrum for r > 0.1 since it does not take into account the fluctuations generated on 
the new last scattering surface. These fluctuations typically do not extend to £ > 100, so if one is 
interested only in scales near the peak of the spectrum a reduction of exp[— 2r] is appropriate. 

Because reionization is such an important part of interpreting CMB anisotropy measurements 
in the context of cosmological models, we will defer a more detailed discussion of reionization to a 



future paper [E^ . 



IV. CURRENT LIMITS ON 1/ 

In this section we discuss what we can infer about u from current observations of CMB 
anisotropy. We shall use the data tabulated in |^ since the newer data is (mostly) not yet available. 
As discussed in and below, the constraints even without including the new data are very 



interesting for large-scale structure modelling. 

Several issues regarding foreground contamination and possible systematic errors intrude in the 
analysis of current CMB data. It was shown in |^ that for the experiments dominating the fit 



near the peak, removing the foregrounds does not increase the error bars by more than ~ 10%. To 
account for foreground removal we have (conservatively) multiplied the errors in |Q by 10% (except 
for COBE). Note that since several experiments have reproduced earlier observations, indicating 
that systematic errors are not as severe as might have been thought. In one case however, that of 



the MSAM experiment |32], one of the two channels (the "single-difference") showed a discrepancy. 



7 



To be very conservative we have dropped this point from our analysis (shown in Fig. ^ as the sohd 
square at £ ~ 150) and updated the "double-difference" point to include the new data (which is 
slightly lower than the older data). Also we do not include the new data from the Python experiment 
(though we show it in Fig. |^ as the solid triangles with the small error bars [^) since the points 
are correlated in an unknown way. 

As it is not the purpose of this paper to revisit the data analysis, especially without access 
to the latest data, we shall use the data set as tabulated in with the simple modifications 
discussed above. One point deserves special mention. In |^ we used symmetric error bars on all 
of the points. This is a conservative method if the inference is the presence of a peak, since most 
of the observations have skew positive error bars. However this leads to a slightly stronger than 
warranted upper limit on the height of the peak, and a slightly weaker lower limit. As most of the 
skewed likelihood functions have not been published in tabular or graphical form we have not tried 
to correct for this bias. Since for most of the models under consideration (especially the models 
with high Q\) the large-scale structure data is best fit with a tilt (n < 1), this treatment remains 
the most conservative. 

A fit to the data as in 0, but with spectra generated from Eq. |^, gives a likelihood function for 
ly shown in Fig. ^. Here we have integrated, or marginalized, over the normalization dIq^, which is 
well fixed by COBE []I0[. The mean and standard deviation are u = —0.05 ± 0.07. The absolute 
goodness of fit to both the COBE data alone |]18[ and the other data [^] show no indications that 



the fit should be suspect for purely statistical reasons. 

Notice that this result has several immediate implications for structure formation models. That 
the preferred peak height is near that of a CDM model with scale invariant initial conditions puts a 
lower limit on the amount of tilt which can be accommodated P,|5(]|. The lower limit clearly depends 
on the values of Qq, h and Qsh"^ assumed, with lower Qq and h acting to loosen the constraint and 
lower flBh'^ acting to strengthen it. The probability of early reionization ||25| in these models will 
only tighten the lower limit on n. 

The increased height of the peak in Iow-Qq models or models with < 1 can cause conflict with 
the upper limit on z/, unless the models are tilted, the baryon fraction is lowered or there is some 
reionization. We note in passing that for ACDM, some tilt is probably necessary to provide a good 
fit with the large scale structure data in any case P, |l8| , ^p5| . In general, while the degree scale 
data provide a strong constraint on models with high Qq, the possibility of early reionization (the 
degree scale power is exponentially sensitive to the very uncertain zr) and the large uncertainty in 
flBh"^ from BBN significantly weaken the constraints from this data for low Qq or models. 



V. CONCLUSIONS 

In this paper we have compared the accuracy of spectra generated using the "cosmic confusion" 
assumption to those calculated using numerical evolution of the Boltzmann hierarchy. For a 
range of parameters of current interest the spectra are found to agree to ~ 10% for all multipoles 
up to the first peak (£ ~ 220). Over the same range the height of the peak is reproduced to ^ 5%. 
As pointed out by Dodelson & Kosowsky [^], the presence of this degeneracy for scales greater than 
0°5 allows a simple statement of the parameter constraints arising from CMB data. Given current 
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data we find that v (see Eq. [TUp is constrained to be = —0.05 ± 0.07. By limiting the amount 
of spectral tilt, this result in combination with large-scale structure measurements (on somewhat 
smaller scales) strongly constrains models of structure formation (see e.g. [^,^,^] and §fV[). 
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TABLES 



I 






Di 


I 


Di 


2 


0.91 


54 


1.63 


178 


4.63 


4 


0.92 


64 


1.78 


195 


4.93 


6 


0.95 


74 


1.94 


212 


5.07 


10 


1.00 


84 


2.12 


230 


5.03 


14 


1.06 


96 


2.38 


248 


4.80 


19 


1.13 


108 


2.67 


267 


4.39 


25 


1.22 


120 


3.00 


287 


3.87 


31 


1.31 


134 


3.41 


308 


3.29 


38 


1.41 


148 


3.83 


329 


2.79 


46 


1.52 


163 


4.26 


351 


2.42 



TABLE I. The multipole moments, Df, = £{e + l)Ce, for the "standard" Cold Dark Matter model with 

(T) 

= 1, 0,B = 0.05, n = 1 and =0. The moments have been normalized to have Diq = 1. Enough 
moments have been given to interpolate the behaviour of the curve. 
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h 


n 


Ub 


zr 


r 


0.3 


0.90 


0.92 


0.02 


28 


0.165 


0.4 


0.70 


0.93 


0.03 


23 


0.141 


0.5 


0.57 


0.94 


0.05 


20 


0.126 


0.6 


0.49 


0.95 


0.06 


18 


0.115 


1.0 


0.45 


0.80 


0.10 


11 


0.065 


1.0 


0.50 


1.00 


0.05 


45 


0.271 



TABLE II. A "reasonable" estimate of the ionization redshift, zr, for some models. The optical depth 
r from z = to zr assuming full ionization of the hydrogen is also shown. The models with < 1 have 
been chosen to be "best fits" to LSS data. Specifically they have = 0.015, a "shape parameter" 

r = 0.25 and as = 0.70^°"^ (which fixes n). Since n < 1 the effective shape measured by the galaxy 
correlation function will be TeS c± 0.21-0.24. The last two models are a "best-fitting" and a "standard" 
CDM model. 
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FIG. 1. The dimensionless amplitude of the matter power spectrum at horizon crossing, 6h, vs. matter 
density = 1 — Cl\ fov n = 0.85,0.9,0.95, 1.00 as fit to the 2-year COBE data. Sh is defined in terms 
of the r.m.s. fluctuation per logarithm in wavenumber, k, by A^(A;) = k^P{k) /{2'it'^) = S'fj{k/H)^~^^T'^{k). 
The solid line is a scale-invariant spectrum: n = 1. The dotted lines show values for lower n (lowest n 
is highest line) with no tensor contribution. The dashed lines show the same range of n (now lowest n 
is lowest line) for models with a tensor contribution as predicted by power-law inflation. These results 
include the dependence of the tensor-scalar ratio on Q,\ that was incorrectly omitted in Ref. 18. 
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FIG. 2. The CMB anisotropy power spectrum, i{e+l)Ce vs log^, for the "standard" Cold Dark Matter 

(T) 

model with = 1, i^s = 0.05, n = 1 and — 0. The vertical axis is approximately the power per 
logarithmic interval in multipole i, which probes angular scales 9 ~ £~^. For experiments on angular scales 
> 0°5 we are interested only in the rise into the first peak (at £ ~ 220). 
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FIG. 3. The CMB anisotropy power spectra for two models with the same Uoh^ = 0.1 and 
rifi/i^ = 0.0125. The sohd Hne is for = 1 and = 0.1, while the dashed line is fio = 1 — = 0.4 and 
h? = 0.25. Notice that the spectrum is not degenerate in Q.ohP' and Q,Bh^ at the 5-10% level, even up to 
the first peak. 
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FIG. 4. The current situation with respect to CMB data, which provide constraints on v. The open 
circles come from a maximum hkelihood fit to the COBE 2-year data (Ref. 35) and arc highly correlated. 
The filled circles are the experiments tabulated in Ref. 5. The errors are ±la, but the upper limit is 
95%CL. The horizontal bars on each point show the range of scales probed by each experiment. (See 
Ref. 5 for more details.) The points which are labelled by squares and triangles are not included in the 
analysis: see text for discussion. 
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FIG. 5. The likelihood function for adapted from Ref. 5. This likelihood has been integrated (or 

1 /2 

marginalized) over the normalization D^q , which is well fixed by COBE. 
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